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SUMMARY 


Reviewed in this paper are studies carried out and/or in progress at the 
MIT Aeroelastic cind Structures Research Laboratory to develop theoretical pro- 
cedures for predicting the large-deflection elastic-plastic transient structural 
responses of metal containment or deflector structures to cope with rotor-burst 
fragment impact attack* Most of the past effort was devoted to containment/ 
deflector (C/D) structures whose axial dimension is comparable to that of the 
attacking fragments and hence the associated structural responses are essen- 
tially two-dimensional* Recent effort has bean applied to analyzing C/D struc- 
tures whose '"axial dimension" is much larger than that of the attacking fragments; 
thus, the associated structural response to be analyzed is essentially three- 
dimensional* 

For two-dimensional C/D structures both finite-element and finite-difference 
analysis methods have been employed to analyze structiural response produced by 
either (a) prescribed transient loads or (b) fragment impact* For the latter 
category, two time-wise step-by-step analysis procedures have been devised to 
predict the structural responses resulting frcxn (a succession of) fragment impacts: 
(1) the collision force method (CFM) whereby one utilizes an approximate predic- 
tion of the force applied to the attacked structure during fragment impact (also 
equal and oppositely to the fragment itself) and (2) the collision imparted veloc- 
ity method (CIVM) in which one c<»^>utes the in^ct- induced velocity increment 
acqiured by a region of the impacted structure near the impact point (and the 
.attendant velocity decrement suffered by the attacking fragment). The merits and 
limitations of these approaches are discussed* For the analysis of 3-d respo ^es 
of C/D structures, only the CIVM approach is being investigated. 

Experimental data for assessing the accuracy, limitations, and versatility 
of these analyses have been obtained from two sources* The Naval Air Propulsion 
Test Center has provided data on the responses of contairu«.ant rings to (a) a 
single T58 turbine rotor blade and (b) to tri-hub burst fragment attack from a 
T58 turbine rotor. Simpler impact experiments involving a "non-deformable 
fragment" (a solid steel sphere) against simple aluminum beams and panels have 
been conducted at the MIT Aeroelastic and Structures Research Laboratory, Com- 
parisons of predictions with observed structural response data are discussed. 
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under N6R 22-009-339; the authors wish to acknowledge also the help of their 
various colleagues (see co-authors cited in references of Appendices A and B), 
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Introduction 


Engine rotor burst fragments may impact against the engine casing and/or 
against special protective structures* These structures ma^ be intended either 
to contain or to divert the fragment and to allow it to escape along a "harmless” 
path; the respective behavior is termed as being either fragment containment or 
fragment deflection. Of principal interest in this paper is the theoretical 
prediction of container or deflector structures (C/D structures) which are 
subjected to fragment intact. Further, attention is restricted to single- 
layer metallic protective structures; the use of non-roetallic materials for 
protective structures tindergoing fragment intact is addressed by several other 
papers in this Workshop* 

If the dimension of the protective structure in the direction parallel to 
the axis of rotation of the turbojet engine is cos^>arable to the corresponding 
dimension of the attacking fragment, the deflection of the attacked structure 
will be essentially the same at all locations along that axial direction; in 
this case, the deformation is termed two-dimensional (2-D). However, if that 
protective-structure dimension is large in the above c<»^;arative sense, the 
structure will undergo general three-dimensional (3-D) structural deflections. 

For preliminary design and parametric studies of C/D structures, it may 
be useful to idealize the transient structural response as 2-D, as depicted 
schematically in Fig. 1. Here the effect of the structure %diich supports the 
C and/or D structure is represented by a normal and tangential spring founda- 
tion; also, various support conditions can be provided in this type of idealized 
2-D model. This type of model tends to include the main structural response 
features while minimizing the con^utational burden. Accordingly, a series of 
2-D structural response codes for partial and/or cooplete rings of arbitrary 
initial shape, with uniform or nonuniform thickness, and subjected to initial- 
velocity distributions, prescribed extemally-appli d loads, or fragment impact 
have been developed. The capabilities and features of these con^uter codes 
[1-5]* are summarized in Appendix A. Some illustrative examples of the use of 
some of these codes are shown later in this paper. 

For structural response conditions wherein the use of a 2-D idealization 
is an excessive over-siiq>lification and %diere one seeks to predict the response 
in greater detail, the structure needs to be modeled as an asse8dt>lage of shell 
elements (and stiffeners) [6-8] to enable an accounting of the 3-D shell 

References are Indicated by numbers in square [ ] brackets. 
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structural deflections which are present. On the other hand an excessively 
fine modeling such as the use of 3-D solid elements to represent a singl'^.- 
layer shelly stiffeners, etc. leads to an excessive computational burden for 
many purposes. Hence, "shell behavior” modeling serves as a logical "next 
improvement" over 2-D modeling of C/D structures. Accordingly, theoretical 
prediction methods to compute the responses of plates and shells to initial 
velocity distributions and prescribed externally-applied transient loads [6] 
are being adapted to predict structural response to fragment impact (91. 

In order to evaluate the accuracy and adequacy of these structural 
response prediction methods, various experiments have been carried out. The 
Naval Air Propulsion Test Center (NAPTC) has provdied data on the responses 
of aluminum and steel containment rings to (1) impact by a single T58 turbine 
rotor blade and (2) to tri-hub burst fragment attack from a T58 turbine rotor 
([10-13], for example); in these cases the attacking fragment is complex and 
undergoes a considercd>le amount of deformation during its impact interaction 
with the containment ring. A cleaner, less-complex set of impact experiments 
has been conducted at the MIT Aeroelastic and Structures Research Laboratory, 
involving steel-sphere intact against (1) beams, (2) uniform-thickness initially- 
flat square aluminum panels, and (3) panels of type (2) but with integral stiff- 
eners of rectangular cross section; transient strain, permanent strain, and 
permanent deflection data of good reliability and accuracy for c<xnparison with 
predictions were obtained [14,15]. Some of these studies are described briefly 
in the following. 

At the present time, theoretical-experimental correlation studies utilizing 
the NAPTC and the MIT-ASRL experimental data are in progress for the 2-D cases; 
for these cases the CIVM-JET 4B computer code is being en^loyed. For fragment- 
inqpact panels (which undergo 3-D responses), same preliminary calculations are 
under way using the breadboard CIVN-PLATE code; systematic testing and checking 
of this code will be required before it can be used with confidence. 

Figure 2 serves as a concise outline of most of the MIT-ASRL studies which 
have been carried out to date concerning the theoretical prediction of the 
responses of metallic C/D structures to fragment inq>act. Listed in Appendix B 
are the associated MIT-ASRL reports and papers, as well as the status and 
availability of the pertinent structural response computer codes. 



Section 2 is devoted to describing two of the analysis methods (the 
collision imparted velocity method CIVM^ and the collision force method CFM) 
studied for predicting the large-deflection, elastic-plastic, transient 
responses of 2-D structures which are subjected either to impulse loading or 
to fragment impact attack; illustrative excunples of the application of these 
methods are shown, with emphasis on the CIVM approach. Section 3 deals with 
theoretical and experimental studies of fragment- impact- induced responses of 
panels which undergo 3-D structural responses. Comments are given in Section 4 
concerning the status of structural response prediction procedures for 2-D and 
3-D single-layer metallic C/D structures, as well as observations concerning 
analysis needs for multilayer multimaterial C/D concepts and configurations 
being considered as lighter weight candidates to cope with energetic engine 
rotor burst fragments. 

2. 2-D Structural Response Studies 

Some representative analyses and results will be illustrated here concisely 
more extensive results and discussion may be found in the cited references. 
Analysis of, 2-D structural response to fragment impact will be discussed for two 
approaches: the collision imparted velocity method (CIVM) and the collision 

force method (CFM); for illustration, both approaches are applied to analyze 
the transient response of a containment ring to impact by a single blade of a 
turbine rotor. Next, a more complex fragment attack is analyzed hy using the 
CIVM approach; this involves T58 tri-hub turbine rotor burst attack against a 
steel containment ring. Because of the complexities arising mainly from severe 
changes In the geometries of the attacking fragments during the impact and 
interaction process, it became advisable to obtain experimental data for a 
more clearly defined impact situation in order that the measured transient 
response information could be used to make a clear assessment of the adequacy 
of the basic building blocks contained in the theoretical prediction procedure. 
Accordingly, described next are experimental and theoretical studies of the 
transient responses of single beams to in^ulse loading or to steel sphere 
intact attack. 

2.1 Single Rotor Blade Impact Against a Containment Ring 

In these studies of 2-D structural response to impact, use is made of 
finite element and finite difference methods %diich have been showit to produce 


reliable predictions for large-deflection, elastic-plastic, transient response 
of simple beams and rings subjected to known impulsive loading [16,17J, For 
impact- induced structaral response analysis, the principal added ingredient 
to be taken into account is the impact/interaction itself; two methods (CIVM 
and CFM) explored for treating this matter are discussed next* 

To illustrate these approaches, impact of a single blade from a T58 turbine 
rotor against a containment ring will be studied* High speed photographic data 
for such a case have been obtained at the spin-chamber facility of the Naval 
Air Propulsion Test Center [10] . Ring configuration data and blade orientation 
as a function of time at intervals about 30 microseconds apart were obtained 
and are used for illustrative comparisons. 

2.1.1 Analysis with the Collision Imparted Velocity Method 

Figure 3 illustrates a containment ring which is modeled by a number of 
finite eluents and subjected to impact attack by an idealized single rotor 
blade. The equations of motion for the ring and for the fragment are solved 
in small increments At in time by an appropriate finite-difference time operator 
sche ^* For this analysis, the fragment is regarded as being rigid. lo^ct is 
regarc i as being an instantaneous local effect between the fragment and a small 
region of the structure in the vicinity of the inq>act point; for present purposes# 
the size of this small ring region on either side of the inq>act point is estimated 
as being the product of At and the longitudinal elastic wave speed in the ring. 
lflq)ulse/momentum and kinetic energy conservation equations are used to calculate 
the “post-inq>act” (or collision- imparted) velocities of the fragment and of this 
in^ct-affected structural region; by employing the concept of the coefficient 
of restitution (e) this local impact can be treated as perfectly elastic (e=l), 
perfectly inelastic (e=0), or intermediate (0<e<l). 

Figure 4 is an information flow diagram illustrating the use of this 
'*collision-iB^arted velocity method” (CIVH) in the calculation of transient 
structural response produced by fragment inq>act* Typically, a succession of 
impacts is predicted* Fuller details of this approach are given, for example, 
in Refs* 4, 17, and 18* 

Table 1 summarizes the containment ring and fragment data for the illustra- 
tive case: NAPTC Test 91* Shown in Fig* 5 are predicted and measured deformed 

ring configurations and blade locations at 150, 570, and 810 microseconds after 
initial impact. In the ii^ct quadrant the ring was modeled by 10 equal length 



cubic-cubic elements; 6 equal-length elements were used in each of the 3 other 
quadrants. For these calculations, coefficients of restitution of e=0 and e=l 
were used. The 6061-T6 aluminum ring material was regarded as elastic, £erfectly- 
£lastic, with a £train-£ate (EL-PP-SR) dependent yield stress given by the 
following approximation: 


1 



where a is the static yield stress, : is the strain rate, and D and p are 

^ -1 
material constants. For these calculations, 0^6500 sec and p=4 were assumed. 

Also, frictionless impact and interaction (y=0) between the ring and the blade 

was assumed for the cases illustrated here* Fairly good agreement between the 

predicted and observed deformed ring configuration is noted, but the predicted 

vs* observed fragment motion is not good* This latter disagre^ent stems mainly 

from ignoring friction and the changing mass moment of inertia of the actual 

deforming blade. Later calculations included these effects. 

2.1.2 Analysis with the Collision Force Method 

In this method the attacking fragment is treated as being deformable. 

Shown, for example, in Pig. 6 are some postulated idealized configurations to 
represent a deformable impacting blade* The straight rigid blade model was 
used in the previous case* Es^lored in Ref* 19 were the following two idealiza- 
tions — the blade was assumed (a) to r^nain straight but to shorten in an 
elastic, perfectly-plastic (EL-PP) fashion or (b) to curl in a simple plausible 
assumed-mode fashion; these are termed, respectively, the elastic, perfectly- 
plastic shortening blade model (EL-PP-SB) and the elastic, perfectly-plastic 
curling blade model (EL-PP-CB) • These modes of behavior combined with a step 
by step collision inspection set of rules permitted following this process. 

At any given instant, applicable values of governing geometric deformed-.* ade- 
configuration parameters were identified. These in turn were related via 
energy methods to the coiq>onent of the force applied by the blade perpendicular 
to the surface of the attacked containment ring and equal-and-oppositely to the 
fragment itself. Similarly, equal and opposite tangential forces (from friction) 
were posttdated to be y times the normal- to- the-sur face component. A self- 
explanatory information flow chart for the CFM process is given as Fig. 7. 

Shown in Fig. 8 are defomed ring predictions at two instants after initial 
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impact for the EL-PP-SB model and the EL-PP-CB model for the case in which the 
friction coefficient p is assumed to be .15 and the perfectly-plastic yield 
stress of the steel--alloy blade is assumed to be - 160,000 psi. Fairly 
good agreement between experiment and these EL-PP-cl model predictions is 
observed* Reference 19 shows the ring response to be rather insensitive to 
(plausible) values of friction coefficient used. The motion of the blade, 
however, is much more sensitive to p -- as Fig. 9 indicates. 

The curling blade model [19] devised by plausible engineering rules and 
approximations appears to represent rather well the behavior and the observed 
deformed configuration of the actual single blade in ffkPTC Test 91. One must 
keep track of the time-varying ge<xnetry of both the deforming blade and the 
deforming containment ring in order to determine when and where the successive 
collisions (i.e., attempted simultaneous occupamcy of scxne regions of space) 
occur. Hence, it is evident that if one were to use this method to ^malylte 
structural response to impact by, for example, a disk-rim fragment with perhaps 
3 to 10 attached blades (each of which will undergo sequential different deforma- 
tions), one would be faced with a substantial book-keeping job to define the 
space occupancy of this complex deforming fragment; the advisability of seeking 
a less complex scheme is clear. Accordingly, subsequent attention has been 
given to the use of greatly- idealized rigid fragments in conjunction with the 
CIVM analysis scheme. 

2.2 CIVM Analysis of Tri-Hub Rotor Burst Attack Against a Containment Ring 

One type of postulated engine rotor fragment attack which has received much 
discussion is that in %diich the rotor bursts into 3 equal segments (termed a 
tri-hub burst) . One fragment of this type is shown schematically in Pig. 10. 

The NAPTC has conducted many tests involving tri-hub burst attack against various 
single- layer and multilayer containment rings. Recently NAPTC Test 201 involving 
tri-hub burst attack of a T58 turbine rotor at 19,859 rpm against a cast 4130 
steel containment ring of 7. 50- in inner radius, 0.625- in thickness, and 1.50- in 
axial length was conducted (13] « Figure 11 shows the post-test deformed-ring 
configuration. High speed photographs showed the severe deformation incurred 
by many of the blades during the ii^pact/interaction process; this is depicted 
schematically in Fig. 10. 

For convenience and geometric simplicity, each such fragment has been 
idealized for use in the CIVM- JET 4B computer code (4] as a rigid circular body 
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of the same mass and mass moment of inertia as the pre- impact fragment, with 
the same CG location, translational velocity, and rotational velocity as the 
actual fragment at postulated release. As indicated in r’ig. 10, one might 
elect to represent the actual fragment by an idealized fragment of "properly 
selected radius ^ examination of this rotor indicates that reasonable 

minimum and maximum values for r^ would be about 2.56 and 4.20 inches, respec- 
tively; the use of these as well as an "intermediate" value of 3.36 inches was 
explored. 

Figure 12 indicates the geometric, test, and modeling data for this case. 

The ring has been modeled by 48 equal-length ring elements. The point of 
initial impact of each of the three fragments is indicated in Fig. 12; element 
numbers and node identification are also given. The uniaxial static stress- 
strain properties of 4130 cast steel were approximated by piecewise linear 
segments with the stress-strain pairs: (a,e) == 80,950 psi, .00279; 105,300 

psi, .0225; and 121,000 psi, .200 via the mechanical sublayer model; strain 
rate effects were approximated by using D = 40.4 sec ^ and p=5. Shown in 
Fig. 13 is the predicted ring configuration at 1000 micioseconds after initial 
impact. The predicted inner surface and outer surface strains at the midelement 
location of elements 1, 4, and 6 are given in Fig. 14; for this calculation, 
fricti<^nless impact y=0 and r^=2.555-in were employed. Figure 15 shows the 
circumferential distributions of inner-surface and outer-s\irface strain at 
2400 microseconds after initial impact. 

The effects of friction for otherwise identical modeling are indicated 
roughly by the Fig. 16 comparison of deformed ring configurations at 1200 
microseconds after initial impact for y==0 and ii=0.3. Similarly, the effects 
of idealized fragment radius are seen in Fig. 16 where deformed ring profiles 
at 1200 microseconds after initial impact are shown for r^=2. 555-in and r^=3. 360- 
in. It is evident that if one chooses an unduly large idealized fragment radius 
r^, this "rigid fragment" will constrain the ring to restrict its bending strain 
contribution so that unrealistically small total strains will be produced at 
the "convex lobes" — compared with that which the actual "effectively-sroaller- 
radius" fragment will produce. 

The use of an idealized fragment of constant radius will clearly make it 
impossible to obtain complete time history agreement between predicted and 
measured inner-surface and/or outer-surface strains. However, the hope is that 
a properly-chosen effective r^ will lead to reasonable predictions vs. experiment 


of maximum strains produced as a function of circumferential location. Further 
calculations and measurements are neede- i»sess the reliability with which 
this can be done. However, at the cost of greater complexity and computational 
expense, one can devise and use a fragment model which more closely simulates 
the behavior of the actual fragment. 

Note, finally, that a comparison between the predicted and observed 
permanently-deformed ring configuration is not shown. This is the case 
because the calculation at At="l microsecond has been carried out only to 2400 
microseconds after initial impact. Whereas peak response occurred near 1200 
microseconds, the ring is still springing back considerably at the 2400 micro- 
second time. A longer calculation would be necessary in order to permit making 
a reasonable estimate of the permanent-deformation configuration. 

2.3 Beam Response to Steel Sphere Impact 

In order to obtain appropriate and detailed 2-D transient structural 
response data under well-defined impact conditions so that a definitive 
evaluation could be made of the adequacy of the approximate collision- 
interaction analysis employed in the CIVM scheme, some simple experiments 
have been conducted at the MIT-ASRL. Beams of 6061-T651 aluminum with nominal 
8- in span, 1.5-in width, and 0.10- in thickness and with both ends ideally 
clamped (see Fig. 18) have each been subjected to midspan impact by a solid 
steel sphere of one- inch diameter [14] • Impact velocities ranged from those 
sufficient to produce small permanent deflection to those needed for threshold 
rupture of the beam. Spanwise-oriented strain gages were applied to both the 
upper and the lower (impacted) surface of the beam at various spanwise locations. 
In each test, transient strain measurements were attempted for 8 of the gages; 
after each test, permanent strain readings were obtained for all surviving 
gages. Also, permanent deflection measurements were made. 

An inspection of each specimen indicates that except near the point of 
impact itself (i.e., where |x|>0.8-in), the beam underwent essentially 2-D 
deflection behavior; pronounced 3-D behavior occiirs near the point of initial 
impact. Hence, the 2-D structural response code (CIVM- JET 4B) may be expected 
to provide valid comparisons for (x|>0.8-in. Accordingly, such calculations 
and comparisons are in progress, and some preliminary results are shown next. 

For the test and specimen identified as CB-18 in Ref* 14, the entire beam 
has been modeled with 43 equal- length cubic-cubic finite eluents. The beam 
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material has been modeled as having either elastic, strain-hardening (EL-SH) 
or EL-SH-SR behavior where the uniaxial static stress-strain curve has been 
approximated by the a,e pairs: a,e= 41,000 psi, .0041; 45,000 psi, .0012; 
and 53,000 psi, .1000. For EL-SH-SR conditions, D=6500 sec~^ and p=4 have 
been assumed. For CB-18 initial steel-sphere impact occurred at a velocity 
of 2974 in/sec; a state of large permanent deflection was produced. 

Shown in Fig. 19 are predicted and measured strains at spanwise stations 
x=1.50 and 1.20-in from the midspan impact point. At these 2-D structural 
response locations, there is fairly reasonable agreement between predicted 
and measured strains. Figiire 20 shows the predicted transient vertical dis- 
placement response at x=1.0-in for both the EL-SH and the EL-SH-SR case. From 
these and longer-duration plots, the estimated respective permanent deflection 
is 0.63 and 0.58- in; the measured value is 0.60- in. While the comparisons shown 
here indicate encouraging agreement* nk.re extensive calculations and compeurisons 
are needed before a firm assessment can b ' made of the adequacy of the procedure 
embodied in the CIVH-JBT 4B computer code '4] . 


3-D Structural Response Studies 


Of concern here are situations in which the fragment-impacted structure 
undergoes pronounced 3-D rather them 2-D deformation. Appropriate methods of 
structural response analysis ^md corresponding well-defined experimental 
tr£U)sient structural response data which will serve to permit making a clear 
evaluation of the adequacy cuid/or accuracy of proposed prediction schemes are 
needed. Some contributions to this process are described here. 

Although structural response analyses for fragment impact against initially- 
curved as well as initially-flat target structures are of interest, it is useful 
to minimize the complexities while checking the adequacy of the basic building 
blocks in the analysis process. Hence, attention has centered on impulse and 
inq^ct experiments and theoretical analysis of initially-flat stnictures. 
Experiments involving steel-sphere impact against (1) narrow-plate (or beam) 
specimens [14] as well as (2) square uniform-thickness panels with four clamped 
edges and (3) panels of type (2) but with integrally-machined stiffeners of 
rectangular cross-section [15] have been conducted. 

Two of the type (2) initially-flat specimens have been subjected to well- 
defined in^ulse loading by the sheet esqplosive loading technique to produce 
large-deflection, elastic-plastic transient structural response data for checking 
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the basic finite-element and transient response prediction aspects — independent 
of impact itself. Also, steel sphere impact tests of this type of panel have 
been conducted. Thus, transient strain, permanent strain, and pexmanent deforma- 
tion data of high quality are available for checking the prediction procedures 
of Refs. 6 and 9? the latter pertains to the breadboard computer codes PLATE 
and CIVM-PLATE which hopefully will enable one to predict 3-D transient large- 
deflection elastic-plastic, structural responses of panels caused by impulse 
and impact, respectively. If future correlation calculations reveal these codes 
to provide relifidjle transient response predictions, these codes will be upgraded 
to a condition convenient for routine use. 

To illustrate the general character of the panel deformations produced for 
this purpose. Fig. 21 shows the permcuient deflection along the centerline of 
specimen CP-2, a 0.062- in thick square initially-flat 8-in by 8- in panel of 
6061-T651 aluminum with all four edges ideally clamped. The sheet explosive 
loading technique was used to impart essentially a uniform initial normal 
velocity of 16,235 in/sec over a 2- in by 2- in region centered at the panel 
center. Strain gages were also applied at various locations on the non-loaded 
side of the* panel; both transient and permanent strain data were recorded. In 
adf^ition, a pattern of lightly scribed grids was applied to a 3-in by 3- in 
region centered at the panel center on the non-loaded surface. Measurements 
of pre-test and post-test spacings of these grid lines enable one to make a 
rough determination of the permanent relative elongation on that surface as 
a function of location from the center of the pcuiel. Some results from these 
determinations are shown in Fig. 22. 

Steel sphere iit^act against a square 8-in initially-flat 6061-T651 aluminum 
panel of 0.063- in thickness with all four sides ideally clamped results in 
permanent deformation conditions wherein severe permanent deformation is con- 
centrated near the point of initial impact itself as Fig. 23 shows for panel 
specimens CP-8 which suffered 1-inch diameter steel-sphere intact at 2435 in/sec. 
Photo-etched grids spaced 0.020- in apart on the non-in^acted surface permitted 
making the permanent relative elongation measurements indicated in Fig. 24; the 
** large strains” are seen to be concentrated near the impact location and decrease 
rapidly with distance from the center of intact. 

Finally, some illustrative preliminary results from applying the bread- 
board CIVM-PLATE code to steel-sphere-inpacted narrow-plate (or beam) specimen 
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CB-18 are presented here; the steel-sphere impact velocity was 2794 in/sec 
initially. For computational thrift, one quarter of the specimen (see Fig. 18) 
was modeled by a 2 by 11 mesh of flat plate elements .aving 6 degrees of free- 
dom per node, with symmetry conditions imposed along x=0 and y=0; this finite 
element mesh is shown in Fig. 25. Initial impact was assumed to occur at 
(x,y)=(0,0) whereas it actually occurred at about .06-in from this location. 
Relative elongation time histories predicted in this calculation along y=0 
at stations x-0.6-in and x=1.2-in are compared with experimental measurements 
in Fig. 26. Figure 27 demonstrates that this 3-D structural response model 
exhibits 3-D deflection predictions — vertical displacements predicted along 
y=0 (the centerline), y=. 375-in, and y=. 75-in as a function of spanwise location 
x are shown at 800 microseconds after initial impact. The anticipated larger 
displacement is seen to occur along y=0, with decreasing displacements (at 
given x-locations) more remote from the center of impact. Finally, Fig. 28 
shows the predicted lateral transient deflection of the center of the plate 
(x,y)=(0,0) and the observed permanent deflection at this location; reasonable 
agreement is evident* 


Summary Comments 


Presented here is an overview of some of the work carried out to develop 
simple methods for predicting the 2-D transient large-deflection elastic-plastic 
structural responses of metal containment or deflection structures subjected to 
impulse loads or fragment intact; many more details may be found in the cited 
references. This 2-0 type of idealization may serve as a good representation 
of certain fragment/structure impact/interaction situations or as a reasonable 
first approximation to other more con^lex cases. This 2-D idealization is 
relatively inexpensive to apply and may be useful for preliminary design, 
parametric studies, materials screening, etc. Stjructural configurations of 
2-D type included in this discussion consist of con^lete rings, partial rings, 
constant or variable thickness, and uniform or arbitrarily-varying initial 
curvature, with various elastic fo\mdation or various local support conditions 
provided. The associated computer codes (see Appendices A and B) are: 


Structure Subjected to Prescribed Transie n t Loads or Initial 
Velocity Distributions 

JET 3: Single-Layer Structures 

JET 5A: Multilayer Bemoulli-Euler Structures 




structure Subjected Only to Fragment Impact 
CIVM-JET 4B: Single-Layer Structures 

CIVM-JET 5B: Multilayer Bernoulli-Euler Structures 

Comparisons between experiment and predictions indicate good theoretical- 
experimental agreement for JET 3 predictions and a very encouraging but 
incomplete assessment for CIVM-JET 4B; further assessment studies are in 
progress. 

For cases in which the impulsively- loaded structure [6] or fragment- 
impacted structure [9] undergoes significant 3-D structural responses, this 
more complex behavior must be modeled accordingly. Excellent theoretical- 
experimental agreement has been demonstrated [6] for finite-element analysis 
of plates and curved shells which undergo large-deflection elastic-plastic 
deformations in response to known severe impulse loading. Shown in this paper 
are encouraging preliminary comparisons between theory and experiment for 
fragment- impacted structures exhibiting 3-B structural response. Appropriate 
high quality experimental data on steel-sphere- impacted narrow beams, square 
uniform thickness panels, and longeron-stiffened initially-flat panels are 
available for near-future theoretical-erperimental correlation studies to 
assess the accuracy and/or adequacy of the proposed prediction procedures. 

These studies are expected to suggest useful prediction modifications and 
improvements. Extensions to include fragment- impacted ini*-.ially-curved 3-D 
structures would comprise a useful logical addition to the prediction caped^ility. 

Although this discussion has pertained to initially-isotropic metallic 
protective structures, many but not all of these analysis features can be 
carried over to the analysis of multilayer multimaterial protective structures 
— such configurations are of potential future interest, as other papers in this 
Workshop indicate. Although such configurations will be much more con^lex and 
difficult to cmalyze, a validated structural response analysis capability would 
be of considerable value for preliminary design, materials screening, parametric 
studies, and to reduce the amount of ad hoc testing which otherwise would be 
required. The development and checking of accurate prediction methods to 
accommodate structural configurations and materials such as those cited at 
this Workshop in the presentations, for example, of Gerstle, Gardner, amd Holms 
will be a difficult and lengthy process but will represent a highly useful 
state-of-the-art advance. This development and validation will require making 
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careful and detailed transient response observations and measurements for 
well-defined targets (geometry, boundary conditions, material mechanical 
and failure properties) and impact conditions. Realistic types of rotor- 
burst fragments should be used in exploratory experiments and in evaluation 
experiments; such experiments are essential to reveal the principal phenomena 
auid to insure that important response features are not overlooked — as might 
be the case if only highly simplifieu impact experiments were to be conducted. 

However, sin^ler better-defined fragments should be used to minimize uncertain- 
ties when obtaining detailed transient response data which are j.ntended to 
serve as a definitive test of the accuracy cmd/or adequacy of the key building 
blocks of the procedures proposed for predicting the "threshold containment 
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levels" of ^structural responses of multilayer multimaterial C/D structures- # 
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TABLE 1 


DATA CHARACTERIZING NAPTC RING TEST 91 


Ring Data 

Test 91 

Outside Diameter (in) 

17.619 

Radial Thickness (in) 

0.152 

Axial Length (in) 

1,506 

Material 

2024-T4 

Elastic Modulus E (psi) 

10^ 

PP Yield Stress (psi) 

50,000 

Fragment Data 

Type 

T-58 Single Blade 

Material 

SEL-15 

Outer Radius (in) 

7.0 

Fragment Centroid from Center of 
Rotation (in) 

4.812 

Fragment Tip Clearance from Ring (in) 

1.658 

Fragment Length ( in ) 

3.5 

Fragment Length from CG to Tip (in) 

2.188 

Fragment Weight (lbs) 

0.084 

Fragment Moment of Inertia about its 
2 

*4 

CG (in lb sec ) 

2.163x10 

Failure Speed (RPM) 

15,644.4 

Fragment Tip Velocity (ips) 

11,467. 

Fragment Centroidal Velocity (ips) 

7,884. 

Fragment Initial Angular Velocity (rad/sec) 

1,638.3 

Fragment Translation KE (in lb) 

6,756. 

Fragment Rotational KE (in lb) 

290.3 
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DEVELOPMENT OF PREDICTION METHODS FOR STRUCTURAL RESPONSE 


• TO PRESCRIBED TRANSIENT LOADS OR INITIAL 
VELOCITIES— JET CODES (FINITE ELEMENT) 

• TO FRAGMENT IMPACT— CIVM- JET CODES 

• ANALYSIS OP 2-0 STRUCTURES 

A SINGLE-LAYER RINGS — »-NEAR COMPIETKW 

▲ multilayer RINGS — »'IN PROGRESS 

• ANALYSIS OF GENERAL PANEL RESPCMSE 

▲ SINGIE LAYER— IN PROGRESS 

▲ MULTILAYER —»> NEXT 

H EXPERIMENTS 

• SMALL SCALE SIMPLIFIED IMPACT TESTS AT MIT TO 
SUPPIEMENT COMPLEX FULL-SCAIE TEST AT THE NAPTC 

▲ OBTAIN DATA TO MAKE IN-DETAIL EVALUATION OF 
ADEQUACY OF PREDICTION METHOD 

▲ IMPACT OF STEEL SPHERE AGAINST 

# BEAMS 

# FLAT SIN6LE-IAYER PANELS 

# MAFPI£-STIFFENED PANELS 


CLAMPED EDGES 


H THEORETICAL-EXPERIMENTAL C0RRELATI(»1 STUDIES 

• USE OF MIT-ASRL EXPERIMENTS 

• USE OF NAPTC DATA 

B COMPUTER COBBS 

• PARAMETRIC AND SCREENING STUDIES 

• TO ASSIST PRELIMINARY DESIGN OF CCWTAINERS 
AND DEFLECTORS 

FIG. 2 SUMMARY OF MIT-ASRL STUDIES ON ENGINE ROTOR FRAGMENT IMPACT 
ON C/D STRUCTURES 
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o EXPERIMENT 

X CASE CR'llB (EL*PP-SR, e • 0) 

A CASE CR>10B (EL-PP-SR, e • 1) 

RING BEFORE INITIAL IMPACT 



(c) TAXI > 810 Bsac 
Fir. 5 OMCLUlffiO 
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ELASTIC PLASTIC 
CURLING 



STRAIGHT DEFORMING 
ELASTIC-PLASTIC 



NON-DEPORMABLE 


FINAL 



FIG. 6 SCHEMATICS OF ACTUAL AND IDEALIZED FRAGMENTS 



















Z-Coordlnate Axis (in.) 


PRE-IMPACT PROFILE 

EXPERIMENT 

0 EL-PP-SB MODEL 
A EL-PP-CB MODEL 
(Tj. = 0.3 IN) 



(a) First Quadrant tf the Deformed 
Ring at TAII » 326 ysec. 


PIG. 8 COMPARISON OP CPM PREDICTIONS POR EL-PP-SB AND EL-PP-CB BLADE 

MODELS POR y».15 BLADE/RING IMPACT VERSUS EXPERIMENTAL DEPORMED 
RING DATA 
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FIG. 10 SCHEMATICS OF PRE-TEST AMD DEFORMED TRI-HUB DISK/BLADE 
FRA(»IENTS, AMD IDEALIZED MODELS 
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FIG. 11 POST-TEST C(MPIGURATION OF THE STEEL CONTAINMENT RING 
SUBJECTED TO T58 TURBINE ROTOR TRI-HUB BURST IN NAPTC 
TEST 201 



TRI'HUB T58 ROTOR BURST IN NAPTC TEST 201 
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FIG. 13 PREDICTED DEFORNED RING CONFIGURATION AT 1000 MICROSECONDS AFTER 
INITIAL IMPACT FOR THE NAPTC TEST 201 CONTAINMENT RING 




46 BUKENTS RING ELEMENT 1 



PIG. 14 PREDICTED TRANSIENT STRAIN ON THE NAPTC TEST 20 \ CONTAINMENT RING 


RING ELEMENT 4 
MXDSPAN STRAIN 



FIG. 14 eONTlNUKD ^NAPTC TEST 201 RING) 



FIG. 14 CONCLUDED (NAPTC '•.EST 201 RING) 
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PIG. 15 CIRCUMFEREtVTIAL DISTRIBUTION OP PREDICTED INNER-SURFACE AND OUTER-SURFACE 
STRAIN AT 2400 MICROSECONDS AFTER INITIAL IMPACT FOR THE NAPTC TEST 201 
CONTAINMENT RING FOR r.=2.555 IN AND y=0 
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48 ELEMENTS T53 TURBINE ROTOR 

4 DOF/NODE TRI-HUB BURST 
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FIG, 16 COMPARISON OF PREDICTED DEFORMED RING CONFIGURATIONS AT 1200 
MICROSECONDS AFTER INITIAL IMPACT FOR U=0 AND U=0.3 WITH 
r^=«2.555 IN FOR THE NAPTC TEST 201 CONTAINMENT RING 
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48 ELEMENTS T58 TURBINE ROTOR 

4 DOF/NODE TRI-HUB BURST 

EL-SH-SR 
M=0 

X r^=2.555 IN 
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PIG. 17 


COMPARISON OP PREDICTED DEPORMED RING CONFIGURATIONS AT 1200 
MICROSECONDS AFTER INITIAL IMPACT FOR TWO DIFFERENT FRAGMENT- 
SIZE MODELINGS AND FRICTIONLESS IMPACT CONDITIONS FOR THE 
NAPTC TEST 201 CONTAINMENT RING 
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KJEDICTIONS (43 ELEM, 4 DOF/N) 


\ 







PIG. 19 CONCLUDED (CB-18) 





o 



PIG. 21 PERMANENT LATERAL DEFLECTION OF IMPULSIVELY-LOADED SQUARE 6061-T651 ALUMINUM PANEL 
MODEL CP-2 WITH ALL FOUR SIDES IDEALLY CLAMPED 




o 



PIG. 22 UPPER SURFACE PERMANENT RELATIVE ELONGATION DATA FOR IMPULSIVELY-LOADED PANEL MODEL CP- 
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PIG. 23 PERMANENT LATERAL DEFLECTION OF STEEL-SPHERE IMPACTED 6061-T651 ALUMINUM PANEL MODEL CP- 























f ^ ' 


. / 

J 





■ 

* . 9 ' 


I / X' 


H'* 


jt 


I 


ir^;‘ - 




*1 





:ON DISTRIBUTION AT 900 MICROSECONDS AFTER INITIAL 
CB-18 




FIG. 28 PREDICTED TRANSIENT AND MEASURED PERMANENT DEFLECTION OF NARROW-PLATE SPECIMEN CB-18 
AT PANEL CENTER (x,y)=(0,0) 



APPENDIX A 


SUMMARY OF THE CAPABILITIES OF MIT-ASRL COMPUTER CODES FOR PREDICTING 
TWO-DIMENSIONAL LARGE-DEFLECTION ELASTIC-PLASTIC TRANSIENT 
RESPONSES OF RING STRUCTURES 


This description is intended to provide for the reader a convenient 
tabular summary of the principal features and capabilities of the two- 
dimensional transient large-deflection elastic-plastic structural response 
ring codes JET 1 (Ref. 1), JET 2 (Ref. 2 ), JET 3A-3D (Ref. 3), CIVM-JET 4B 
(Ref. 4), and JET 5A and CIVM-JET 5B (Ref. 5) developed under KaSA NGR 22-009- 
339. 
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The JET 1 code of Ref. 1 pertains to single-layer complete, uniform- 
thickness, initially-circular rings of either temperature-independent or 
temperature dependent material properties. These rings may be subjected to 
prescribed : (a) initial velocities, (b) transient mechanical loading, and/or 

(c) steady nonuniform temperatures. The finite-difference method employed 
in this code had been shown previously (Ref. 6) to provide reliable predic- 
tions for the case of temperature- independent material properties. 

The JET 2 code was written in order to extend this finite-difference 
analysis capability to treat multilayer rings — cases anticipated to be of 
future concern. In the interests of efficiency and the minimization of 
computer storage requirements, temperature-dependent material properties 
and thermal loading features were omitted from JET 2; if these omitted 
features should turn out to be needed urgently, they could be added later* 

Since the JET 1 and JET 2 codes pertained to initially-circular, 
complete rings of uniform thickness whereas there was interest also in 
variable- thickness, arbitrarily curved, partial as well as conqplete rings, 
the JET 3 series codes was developed. To accommodate these latter features 
as well as a variety of types of (1) boundary conditions, (2) elastic- 
foundation supports, and (3) point elastic supports, the more versatile 
finite-element analysis procedure was developed and employed. For efficiency 
and user convenience, four versions of the JET 3 program were developed; each 
version accommodates both complete rings and partial rings. JET 3A and JET 3B 
pertain to. uniform-thickness, initially-circular rings, and employ, respectively, 
the central-difference and the Hoxibolt finite-difference time operator; for 
certain cases, the latter finite-difference time operator may permit more 
economic converged transient response predictions than the former. The codes 
JET 3C and JET 3D are corresponding codes which acconmodate variable-thickness, 
arbitrarily-curved rings. 

In most of these codes (JET 1 through JET 3D and JET 5A), the stimuli: 

(1) initial velocity or in^ulse conditions and/or (2) transient mechanical 
loading must be prescribed by the user or analyst. The externally-applied 
forces experienced by a complete or a partial ring from fragment impact are 
not provided within these codes. The user must supply his own estimate of the 
distribution and time histories of these forces. However, in the CIVM-JET 4B 
and CIVM-JET SB codes, fragment/ring interaction and response effects are 
handled internally automatically, for the idealized single-fragment and 
n-fragment cases provided and discussed in the Appendices of Refs. 4, 5, 
and 7. 
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The CIVM-JET 4B code (Ref. 4) was developed from a modified version 
of the JET 3C codCf using the central difference timewise operator. The 
CIVM (collision imparted velocity method) handles a fragment- structure 
impact as a series of quasi-static momentum transfers between the attack- 
ing fragment and the local-impact-affected portion of the impacted structure. 
The solution proceeds as though a series of impulses has been applied to the 
impacted region of the structure. This code provides strain output at each 
Gaussian station, nodal location, and designated additional points for user 
convenience, and calculates the reaction forces at each constrained degree 
of freedom. Another feature of this code is the ability to accommodate 
branches which are used as additional structural supports. These branches 
can have material properties either the same or different from those present 
in the main structure. 

The JET 5A and CIVM-JET 53 codes (Ref. 5) were written in order to 
extend the capabilities of the JET 3D and CIVM-JET 4B codes to multilayer 
structures which are assumed to be hard-bonded and to deform in the Bemoulli- 
Euler fashion. Both codes contain the Houbolt timewise operator and all the 
additional strain emd reaction force output and structural support capabilities 
utilized in the CIVM-JET 4B code. 

In convenient tad>ular form, the principal features and capcdiilities of 
the codes JET 1, JET 2, JET 3A-D, CIVM-JET 4B, JET 5A, and CIVM-JET 5B are 
given in the following: 
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STRUCTURAL RESPONSE COMPUTER CODE STATUS 


CIVH-JET 4B 


Capability Status 

2-D Single -Layer Beams and Complete 

Rings Subjected to Prescribed (Ref. 8) 

Transient Loads or Initial 
Velocity Distributions (No 
Fragment Impact) 

2-D Single-Layer Beams and Complete 

Rings Subjected Only to (Ref. 14) 

Fragment Impact 


Availabilit\ 


2-D Multilayer Bemoulli-Euler 
Beams and Rings Subjected to 
Prescribed Transient Loads or 
Initial Velocity Distributions 


Complete 
(Ref. 15) 


CIVH-JET SB 2-D Multilayer B-E Beams and 

Rings Subjected only to 
Fragment Impact 

PLATE and 3-D Single Layer Initially- 

ClVM-PLATE Flat Panels Subjected, Respec- 

tively, to (1) Prescribed 
Transient Loads and/or Initial 
Velocity Distributions or (2) 
Fragment In^ct Only 


Ccxnplete 
(Ref. 15) 


In Progress 


a: Available from COSMIC, Barrow Hall, University of Georgia, 

Athens, GA. 30601; contact MIT for errata. 

b: Available under a copyright licensing agreement from NIT. 

Contact Prof. E.A. Hitmer, Room 41-219, MIT, Cambridge, 
Hass. 02139. 




DISCUSSION 


B>L, Koff, GE-Cincinnati 

We have conducted tests ard managed to collect pieces of blades that 
were deliberately failed to undei , tand containment ring behavior. It is 
quite obvious that you don't get a three-lobe shape in the ring, because 
as soon as the ring starts deforming locally, all of the other blades in the 
rotor act as a bearing for the ring. This tends to keep the ring round, not 
three-corner or some other shape, by adding quite a bit of support to the ring. 
It suggests that there is more to be lean^ed from the tests you are now running 
on panels, than in oversimplified tests run with a ring that is not supported 
in a manner similar to the engine. When you start adding other support, you 
might find that these simplified panel tests, in fact, more nearly duplicate 
what actually happens, than an oversimplified test with rotor burst fragments. 

E.A. Witmer, MIT-ASRL 

I think it would be very useful in this whole program, if we could have 
people like you, who could suggest to \xs proper models to use for supported 
structures, so that we simulate things in the right way. It's an excellent 
idea. 


As I understand your described tests# you released 1 or 2 blade portions 
from a rotating fully-bladed rotor to intact a containment ring. Similar tests 
done at the NAPTC show behavior very similar to what you describe; the initial 
impact causes the ring to deform and then it comes in contact with the blades 
still attached to the spinning rotor. These blades also deform but do "support” 
the ring and tend to restrain it from deforming as severely as it would if a 
"free ring" were impacted only by the initial attacking fragments. 

A. Weaver, P&W 

As I understand this model, it does a fairly representative job of 
modelling deflections in simple structures, whether they are panels or rings. 
However, it doesn't get at the meat of the containment problem as I see it, 
which is failure. I don't always care about deflections, but I do care when 
and where the ring is going to fail, and how to model that. The 2-D analysis 
con^letely ignored the localized effects going on at the center of in^ct, 
which I believe are very important. 

E.A. Witmer, MIT-ASRL 

You're perfectly correct, there are 3-D effects present where failure 
initiates in the cited beam experiments, and 2-D is clearly an idealization. 

It's a convenient scheme to us to obtain some crude estimates but it certainly 
doesn't address the real problem. The 3-D problem is the important one. For 
the beams and rings discussed here, the structural response behavior is of the 
2-D type essentially everyidiere on the (narrow) rings and also everywhere on 
the steel-sphere- injected beam specimens except near the "impact point" itself 
where 3-D effects are very prominent. Here at threshold rupture, a multiaxial 
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strain state involving very large strains exists. For such regions/ the 
analysis must accommodate large strain plasticity effects and an appropriate 
"failure strain" or similar criterion. This is a matter that is receiving much 
attention now by various groups. 

D. Oplinger, Army-AMMRC 

Is it realistic to assume that you're going to get a structural problem 
rather than a penetration problem? Some of the velocities I saw were fairly 
low; they were a couple hundred feet per second, but when you get up to a 
thousand feet per second, you've got to treat the penetration problem first 
and then you can treat it as a structural response problem. 

E. A. Witmer, MIT-ASRL 

As I understood the fragment velocities cited, they represent the frag- 
ment tip and/or the CG velocities; not the velocity component perpendicular to 
the impacted surface at impact. For rotors with typical small clearance, the 
typical impact angle is very shallow — somewhere in the vicinity of 20-25 
degrees. Hence, for many cases, the typical normal-to-the-surface velocity 
component at impact might range up to perhaps about 420 fps. Depending upon 
the material properties of the structure being impacted, the subsequent 
behavior could involve "penetration" followed by structural response or could 
involve principally only structural response* For roost of the contair^ient 
structure materials being considered, I believe that the latter is the more 
prevalent case. 

D. Oplinger, Army-AMMRC 

I am not familiar with blade materials but what little I know would lead 
me to believe that it would be unusual to get such large curling as you were 
showing. Is that typical of common blade materials, that they can bend over 
like that without snapping into small pieces? 

E. A. Witmer, MIT-ASRL 

For the small T58 turbine rotor used in many of the NAPTC tests, this was 
the observed behavior. However, for the rotors of the newer larger engines, 1 
will ask Mr. Kbff of GE to respond — he can give a better answer. 

B.L. Koff, GE-Cincinnati 

Some of the blades axe high aspect ratio turbine blades# and are more 
typical of aft end turbine stages. The first stage of the HP has blades of 
low aspect ratio and the first stage of an air-cooled turbine consists of a 
hollow structure which usually fragments into many pieces upon impact. Titanium 
fan blades don't curl very much but break up into pieces. 

S. Sattar, P&W 

I want to remark on the basic philosophy or approach to fragment contain- 
ment design. Would it make more sense for us to step back and ask ourselves 
that if you go through this analysis and you have to determine when these 
con^uter programs will predict penetration, you would have to calibrate them 
against tests? Might it not be easier to take a sinqpler approach to predict 
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whether the fragment will be contained or not? It is a case of strength of 
materials or solid mechanics approach, calibrated against spin-pit or specimen 
tests versus these codes to predict the deflections and strains, and then 
finding out at what strain value will the penetration finally take place — 
which you will calibrate anyway, against some tests. I would like a comment 
on that. 

E.A. Witmer, MIT-ASRL 

Your point is, a valid one, however, I think that if one can afford to 
run experiments on every kind of configuration, material, and so forth, to 
obtain the data you seek, that*s one way of proceeding. There is some hope 
that one need not go that far, but instead one can rely upon more basic material 
property information and methods of structural dynamic analysis (at least for 
simple cases) and have a reasonable prospect of predicting analytically when 
these containment-structure failures should occur. I believe that the 3-D 
structural response studies in progress represent a useful step in that direc- 
tion. 


Now, one can immediately dream up a new case which is too complicated for 
any available analysis to handle properly. In such cases one would have to 
appeal to selected experiments; it seems to me any good organization would 
always do that. 

J.W. Leech, ERDA 

Would you comment on why an aluminum alloy was used for the beam model, 
the panel models, emd the containment ring which was subjected to single-blade 
impact, 

E.A. Witmer, MIT-ASRL 

We used 6061-T6 and 6061-T651 aluminum for these specimens for fabrication 
convenience and because their stress-strain properties cure well known; very little 
strain hardening is present. We approximated these properties by piecewise linear 
segments and used them (via the mechcuiical sublayer model) in the transient 
response calculations. 

Incidentally, the NAPTC had static stress-strain tests conducted on the 
4130 cast steel used in their containment ring tests. As perhaps you noticed, 
we did not show any comparisons between our calculations and the experiment 
for NAPTC Test 201 (T58 tri-hub burst against the steel containment ring) 
because we have not concluded that work. You can see immediately that the 
idealization that we used for the fragment, will give tis no hope whatever of 
predicting in detail the tr£msient response. The hope is that .a realistic 
selection of the idealized (rigid circular) fragment may enadDle us to predict 
the peak response reasonably well, but the actual physical situation is just 
so much different from the idealized model that the fine transient response 
details actually present can not be reproduced by this model. But that's 
really expecting too much of that simple model. Of course, the model can be 
refined. One can devise a more con^licated fragment model — one can put in 
the various curling blades (attached to the disk segment) and let them go ahead 
and curl and follow them; a tremendous amount of bookkeeping would be involved. 
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Hence, we elected to tiy to see the potential of this single rigid-circular- 
fraginent idealization. Also, one could modify this simple circular fragment 
to permit deformations approximating roughly the behavior of the blade/disk 
fragment itself to achieve a still simple, but better simulation of the actual 
attacking fragment. 



